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ABSTRACT 

A substantial fraction of the light emitted from young or star-forming galaxies at ultraviolet to 
near-infrared wavelengths comes from the ionized interstellar medium in the form of emission lines 
and a nebular continuum. At high redshifts, star formation rates are on average higher and stellar 
populations younger than in the local Universe. Both of these effects act to boost the impact of nebular 
emission on the overall spectrum of galaxies. Even so, the broadband fluxes and colours of high- 
redshift galaxies are routinely analyzed under the assumption that the light observed originates directly 
from stars. Here, we assess the impact of nebular emission on broadband fluxes in Johnson/Cousins 
BVRIJHK, Sloan Digital Sky Survey griz and Spitzer IRAC/MIPS filters as a function of observed 
rcdshift (up to z = 15) for galaxies with different star formation histories. We find that nebular 
emission may account for a non- negligible fraction of the light received from high-redshift galaxies. 
The ages and masses inferred for such objects through the use of spectral evolutionary models that 
omit the nebular contribution are therefore likely to contain systematic errors. We argue that a 
careful treatment of the nebular component will be essential for the interpretation of the rest-frame 
ultraviolet-to-infrared properties of the first galaxies formed, like the ones expected to be detected 
with the James Webb Space Telescope. 

Subject headings: Galaxies: high-redshift - galaxies: ISM - techniques: photometric - radiation mech- 
anisms: general 



1. INTRODUCTION 

In young galaxies, or galaxies with active star for- 
mation, massive stars photoionize the surrounding 
interstellar medium, thereby adding emission lines and 
nebular continuum flux to the integrated spectral energy 
distribution (SED). Since this nebular component can 
make up a substantial fraction of the overall light 
observed at rest-frame ultraviolet to near-infrared 
wavelengths, modellers have for a long time stressed 
the importance of considering the contribution from 
ionized gas to the SED when attempting to interpret 
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tribution is expected to grow with the redshift of 
the galaxies studied, since stellar populations are on 
average younger and the star formation activity more 
pronounced at high redshifts than in the local Universe. 
Even so, observers often attempt to derive ages and 
masses of high-redshift objects from optical to near- 
infrared broadband photometry data using spectral 
evolutionary models that only predict the SED of the 
stars themselves. Recent examples incl ude studies of 
star-f o rming or young galaxies at z ~ 2 (|Shaplev et al.l 
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of this approach by estimating the fraction of nebular- 
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to-stellar light emitted for galaxies with redshifts in the 
range z = 0-15, under various assumptions about the 
prior star formation history. 

2. MODELLING HIGH-REDSHIFT GALAXIES 

We model t he S EDs of galaxies using the 
IZackrisson et al.l (|200ll) spectral synthesis code, which 
treats nebular emission by feeding the integrated stellar 
component spectrum for different ages into the pho- 
toion ization model Cloudy, version 90.05 (jFerland et all 
1998). For simplicity, we assume that the nebula is 
spherical and radiation bounded, has a covering factor of 
unity, and has a gaseous metallicity Z — 0.020 identical 
to that of the stars. 

Three different types of model galaxies are considered 
- early-type, late-type and starburst galaxies - differ- 
entiated by their assumed ages and star formation his- 
tories. For the early- and late- type objects, we adopt 
the functional form suggested bv lNagamine et ail (|2006l ) 
for the temporal evolution of the star formation rate: 
SFR(t) oc (t/r) exp(-t/r). While Nagamine et al. find 
that r = 1.5 Gyr and r = 4.5 Gyr give reasonable 
fits to the blue and red peaks in the colour distribution 
of low-redshift galaxies in the Sloan Digital Sky Survey 
(SDSS), we here adopt r = 0.5 Gyr and r = 4.5 Gyr 
for our early and late-type models, i.e. a somewhat less 
prolonged star formation history for early- type galaxies. 
This pushes the bulk of star formation (90% of the stellar 
mass formed) to z > 3 (as opposed to z > 1 for a r = 1.5 
Gyr model) for the early-type galaxy, making it more 
representative for some of the "red and dead" galaxies 
recently discovered at z = 1.5-3 fe.glLabbe et al.l 1 2005 ; 
Kriek et all 120061: iMcGrath. Stockton, k. Canalize] 12007 : 
Stockton et al.ll2008h . " 



For starburst galaxies, we assume a fixed age of 50 Myr 
at all z bs and a constant prior SFR(i). The SEDs of the 
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Fig. 1. — The ratio of nebular to stellar light as a function of age 
in filters B (solid), R (dashed) and K (dash-dotted) for the early- 
(thin lines) and late-type (thick lines) model galaxies at 2 obs = 0. 



early- and late-type model galaxies are evaluated at the 
age of the Universe at the observed redshifts we consider 
here: z b s = 0-15, in steps of Az Q b s = 0.2. Redshifts are 
converted into cosmic ages assuming a ACDM cosmol- 
ogy with tt A = 0.74, tt M = 0.26 and H = 72 km s" 1 
Mpc -1 , in accordan ce with the 3 year WMAP results 
(jSpergel et al.ll2007t ). 

For every z b s , the SEDs of these model galaxies are 
convolved with the filter profiles of Johnson/ Cousins 
BVRIJHK, SDSS griz, Spitzer IRAC at 3.6, 4.5, 5.8, 
8.0 /im and Spitzer MIPS at 24 /itm. For each filter we 
calculate / ne b//stars, the ratio of the net flux from the 
nebular component (the flux emitted by the nebula mi- 
nus the absorbed stellar flux, the latter being essentially 
zero at wavelengths higher than 912 A) to that from the 
stars themselves, as a function of z bs- Since we here 
assume the nebular emission to be caused by processed 
far-UV starlight, flux conservation requires that the the 
bolometric /neb //stare wu l always be lower than unity. 
However, when integrating over broadband filter pass- 
bands at wavelengths longward of the redshifted 912 A 
Lyman limit, / ne b//stars > 1 is possible and does occur. 

For all model galaxies, we adopt a constant metallic- 
ity of Z — 0.020 and a Salpeter initial mass function 
(AN /AM oc Af -2 - 35 throughout the mass range 0.08-120 
Mq). We stress that this model only predicts the SED 
of the stellar component and its associated photoionized 
nebula, and that neither dust extinction nor absorption 
in the gaseous medium along the line of sight are taken 
into account. The impact that such effects may have on 
the predicted / ne b//stars ratios is discussed in Section 4. 

3. RESULTS 

To illustrate the time dependence of the nebular contri- 
bution, we plot /neb //stars in filters BRK for our early- 
and late- type model galaxies as a function of age in Fig.[TJ 
The redshift is assumed to be z b s = 0. Nebular emis- 
sion is seen to dominate over the stellar component (i.e. 
/ncb//stars > 1) at very low ages (~ 10 6 -10 7 yr), and 
to remain non-negligible (/ ne b//stare > 0.1) up to ages of 
~ 10 8 -10 9 yrs, depending on the filter and star formation 
history. When / n eb//stars is plotted against time in this 
fashion, without considering redshift effects, / ne b//stars 



decreases monotonically in each filter as the stellar popu- 
lations grow older. The situation becomes far more com- 
plicated if /neb//stars is instead plotted against the ob- 
served redshift of these model galaxies. 

In Fig. El we plot / nc b//stars as a function of z obs for 
the three model galaxies and all the 16 filters considered. 
Contrary to Fig. [H / ne b //stars here displays a number of 
maxima, most clearly seen in the case of the starburst 
model (blue line). 

What is the origin of these peaks in the / ne b//stars 
ratio? For filters with central wavelengths in the range 
from B to H, two sets of peaks at low and high z G b s 
are clearly visible (e.g. at z Q b s ~ 0-3 and z Q b s ~ 8-11 
in the J band), with a valley of low / nc b//stars values 
inbetween. The wide peak seen at high z b s in filters 
with central wavelengths in the range B-H is caused by 
the Lya emission line. The peaks at low z Q b s come from a 
number of different spectral features. In filters B and g, 
these peaks are mainly due to the H/3 emission line and 
nebular Balmer continuum. These emission processes are 
accompanied by the [O III]AA4959,5007 emission lines in 
filters VRI and riz, and by the Ha line in RI and riz. 
The /nob//stars peaks due to Balmer continuum, H/3, [O 
III] and Ha become broader as they redshift into the 
wider filters at longer wavelengths. The most prominent 
of the bumps seen at z b s > 2 for filters with central 
wavelengths in the range from if to 8.0 /xm is due to Ha. 

For each filter, /ncb/jstars is only plotted up to the 
redshift where the Lyman limit at 912 A enters the fil- 
ter. When this happens, absorption within the nebula 
begins to dominate over emission processes, and galax- 
ies become increasingly more difficult to detect because 
of the very low stellar flux emerging shortward of this 
limit. For computational reasons, absorption shortward 
of the Lyman limit is here treated as a negative flux, 
which means that / ne b becomes negative and drops be- 
low the plotted range at this redshift. This prevents the 
/nob//stars curves to extend all the way up to z bs = 15 
for filters in the wavelength range from B to J. 

The rise in / ne b/ /stars seen at low z Q b s in filters at 4.5- 
24 /im is partly due to a number of strong emission lines, 
but in the 8 and 24 /im filters, there is also a considerable 
contribution from nebular continuum at z b s ~ 0. Even 
though the flux from the nebula can become much higher 
than the flux of the stars at the latter two wavelengths, 
this is in fact of little relevance for observations, since the 
overall SED is likely to be dominated by dust emission 
in these filters at low z bs, and not by stars or nebular 
emission. 

Even if we disregard the unobservable behaviour of 
/neb//stars at 8-24 fj,m. at z b s ~ 0, certain combinations 
of filters and redshifts produce / ne b//stars of up to w 60% 
for all the model galaxies considered. This means that 
estimates of the stellar population mass, based on the in- 
correct assumption that all of the observed flux is directly 
produced by stars, may be overestimated by this fraction. 
Moreover, / ne b //stars changes dramatically from one fil- 
ter to the next for fixed z D b s , which implies that nebular 
emission will alter the colours of the overall SED. As an 
example, consider the difference in / ne b//stars between 
the J and 4.5 /im filters at z D b s = 6 for the starburst 
model. Nebular emission will change this colour by rj 0.4 
mag, which means that the best-fitting age is likely to be 
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Fig. 2. — The ratio of nebular to stellar light as a function of z \^ s for 16 different broadband filters. The lines correspond to the different 
model galaxies considered: starburst galaxy (blue), late- type galaxy (green), and early- type galaxy (red). The various peaks in /neb //stars 
are mainly due to strong emission lines which shift in and out of the filters as a function of redshift. High values of /neb //stars signify 
redshift intervals where treating the observed broadband fluxes as direct starlight may be a bad approximation (see main text for details). 



substantially affected. While a first-order correction for 
nebular emission can be applied to observed broadband 
fluxes by assessing the likely equivalent widths of the 
most prominent emission lines (e.g. Lya, H/3, [O III], 
Ha), high-accuracy interpretations of broadband fluxes 
will require a consideration of the entire nebular SED, at 
least in the case of starburst-type galaxies. 

4. DISCUSSION 

Our estimates indicate that nebular emission may be 
an important, and hitherto largely neglected ingredient 
in models used to extract ages and masses from broad- 
band fluxes of high-redshift galaxies. We note that its rel- 
ative contribution is comparable to that of recent updates 
of the evolution o f thermally pulsating asymptotic giant 
branch stars (e.g. iMarastoa 120051 : iMaraston et alJl2006t 
lBruzualll2007f ) , which have recently attracted a lot of at- 
tention in the high-redshift community. A careful treat- 
ment of the nebular component is therefore likely to be 
very important for the interpretation of the ultraviolet- 
to-infrared properties of the first galaxies formed, like 
the ones expected to be detected with the James Webb 
Space Telescope 5 . 

However, the / ne b//stars ratios presented here repre- 
sent no more than tentative estimates, since there are 

5 http://www.jwst.nasa.gov/ 



a number of mechanisms that could boost or diminish 
the effects of nebular emission, and which are difficult to 
control at the current time. 

Owing to the likely breakdown of the closed-box ap- 
proximation at high redshifts, we have here adopted a 
fixed stellar and gaseous metallicity of Z = 0.020 instead 
of implementing an age-dependent metallicity. While ap- 
propriate for many types of galaxies at z bs ~ 0, the typ- 
ical metallicity must realistically be much lower at earlier 
epochs, even though the exact redshift evolution is diffi- 
cult to assess. The choice of Z = 0.020 ensures that our 
/neb//stars ratios are conservative, since lower metallic- 
ities would imply higher Lyman continuum fluxes from 
the stars and less efficient metal cooling of the gas, result- 
ing in boosted / ne b//stars ratios. This is demonstrated 
in Fig. [3l where we have plotted the redshift evolution of 
/neb//stars in the K band for our starburst model galaxy 
at various lower metallicities. Lowering Z increases the 
relative nebular contribution at esscntally all z h s , even 
giving / n eb//stars > 1 in three different redshift inter- 
vals. The emission lines responsible for the most promi- 
nent /neb//stars peaks in this passband have been indi- 
cated with labels. Since / ne b and /stars respond to metal- 
licity changes in different ways, the exact behaviour of 
/neb//stars as a function of Z and z bs actually turns 
out to be quite complex. As an example, consider the 
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Fig. 3. — The ratio of nebular to stellar light as a function of z b s 
in K for our starburst galaxy model at Z = 0.001 (thick solid), 
0.004 (thin solid) 0.008 (thick dashed) and 0.020 (thin dashed). 
The latter metallicity corresponds to that plotted in Fig. [2] The 
emission lines responsible for the most prominent features have 
been indicated by labels. 



Ha-dominated peak at z bs ~ 2.2, which shows very lit- 
tle metallicity dependence in the range from Z = 0.020 
to Z = 0.004, but increases rapidly from Z = 0.004 to 
Z = 0.001. The Ha luminosity changes relatively little 
between these metallicities, and the sudden increase in 
/neb//stars &t Z = 0.001 is instead caused by a drop in 
/stars related to properties of the adopted stellar evolu- 
tionary tracks. By contrast, the rapid rise of / ne b//stars 
at z obs ~ 3.2 between Z = 0.020 and Z = 0.008, fol- 
lowed by a less dramatic growth at lower Z, is caused by 
a combination of the nonmonotonic metallicity depen- 
den ce of the [O III1/H/? emission-line ratio (see Fig. 17 
in iNagao. Maiolino. fe Marconil l2006h and the drop in 
/ stars at Z = 0.001. 

We have also adopted the conservative assumption of 
a Salpeter initial mass function throughout the redshift 
range z Q b s = 0-15, whereas theoretical arguments sug- 
gest that the initial mass function should turn top-heavy 
at high rcdshifts, possibly introducing very massive Pop- 
ulatio n III stars (for a review, see iBromm fc Larsonl 
l2004h . Because of the high Lyman continuum fluxes of 
such stars, this would increase / ne b//stars at high z b s - 

Our estimates of the / ne b//stars ratios are based on 
the assumption that gas is ionized by stars only. The 
presence of active galactic nuclei would boost the flux 
of ionizing photons, and shock heating of the gas may 
also become important in galaxies undergoing mergers. 
Such effects are expected to increase / no b / /stars in certain 
high-redshift objects. 

Our model assumes that all Lyman continuum photons 
are absorbed by gas close to the galaxy and therefore con- 
tribute to the nebular SED, but this may be unrealistic, 
since dust extinction shortward of the Lyman limit would 
subst antially reduce the hydrogen-ionizing flux (e.g llnouel 
l2001h . and hence decrease / ne b/ /stars- The exact level of 
Lyman continuum extinction is difficult to assess even 



in the local Universe, and even more so at higher red- 
shifts. This treatment also implicitly assumes that no Ly- 
man continuum photons can escape into the intergalactic 
medium, either through holes in the gas (i.e. a covering 
factor below unity) or through a matter-bounded nebula 
(i.e. insufficient gas supply to form a complete Stromgren 
sphere) . This approximation must of course break down 
if galaxies are to contribute to cosmic reionization, and 
a few detections of Lyman continuum leakage have in- 
deed been claimed at both low and high redshifts (e.g . 
Steidel. Pettini fc Adelbergedl200H iBergvall et al.ll2006l: 



Shaplev et al.ll2006l ). Any mechanism that allows Lyman 



continuum photons to escape would lead to a nebular flux 
lower than that predicted here. 

Since Lya photons are resonantly scattered by neu- 
tral hydrogen both inside the source galaxies and in the 
intergalactic medium, the strength of the Lya line be- 
comes very sensitive to the assumed geometry and kine- 
matical properties of neutral hydrogen along the line 
of sight. Since such effects are not taken into account 
here, the Lya contribution to / ne b//stars that is rele- 
vant for observers is likely to be lower than that pre- 
dicted by us. This affects the height of the peak at 
high z bs in the / ne b//stars evolution for filters with 
central wavelengths in the range B-H in Fig. [2] At 
the current time, estimates of the Lya escape fraction 
/esc, Lya vary greatl y. A Lya escape fract ion as low as 
/esc, Lya ~ 0.02 (e.g. iLe Delliou et al.ll2005l ) would make 
the impact of Lya on / ne b//stars completely insignifi- 
cant, whereas L n .i,vn ~ 0.8 (e.g. iGawiser et afl 120061 : 
iKobavashi et al.ll2007 ) would imply just a slightly smaller 
effect than that indicated in Fig. [2] Either way, the in- 
creased Lya opacity at epochs p rior to complete cosmic 
reion ization (z Q b s > 6; see e.g. iFan. Carilli fc Keating! 
20061 for a review) should make the Lya flux smaller 
than assumed here. Hence, the / ne b//stars peaks seen at 
rcdshifts higher than this in panels J and H in Fig. [2] are 
likely to be too prominent. 

Whereas a uniform, foreground, scrcenlike distribu- 
tion of dust would decrease the nebular and stellar 
fluxes by an equal factor, thereby leaving / ne b//stars 
unchanged, selective extinction, i.e. the fact that the 
young stellar population and the associated ionized gas 
may have different dust properties than the older stars, 
could in principle drive / ne b//stars either way. Ob- 
servations of starburst galaxies suggest that the nebu- 
lar compo nent may suffer more extinction than the ol d 
stars (e.g. ICalzetti. Kinnev fc Storchi-Bergmannlll994l ). 
but the opposite effect has also been advocated to ex- 
plai n the equivalent widths of Lya in high-rcdshift galax- 
ies (jFinkelstein et al.ll2007h . 
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